Search for: All records

Abstract Climate change is rapidly altering composition, structure, and functioning of the boreal biome, across North America often broadly categorized into ecoregions. The resulting complex changes in different ecoregions present a challenge for efforts to accurately simulate carbon dioxide (CO₂) and energy exchanges between boreal forests and the atmosphere with terrestrial ecosystem models (TEMs). Eddy covariance measurements provide valuable information for evaluating the performance of TEMs and guiding their development. Here, we compiled a boreal forest model benchmarking dataset for North America by harmonizing eddy covariance and supporting measurements from eight black spruce (Picea mariana)-dominated, mature forest stands. The eight forest stands, located in six boreal ecoregions of North America, differ in stand characteristics, disturbance history, climate, permafrost conditions and soil properties. By compiling various data streams, the benchmarking dataset comprises data to parameterize, force, and evaluate TEMs. Specifically, it includes half-hourly, gap-filled meteorological forcing data, ancillary data essential for model parameterization, and half-hourly, gap-filled or partitioned component flux data on CO₂(net ecosystem production, gross primary production [GPP], and ecosystem respiration [ER]) and energy (latent [LE] and sensible heat [H]) and their daily aggregates screened based on half-hourly gap-filling quality criteria. We present a case study with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to: (1) demonstrate the utility of our dataset to benchmark TEMs and (2) provide guidance for model development and refinement. Model skill was evaluated using several statistical metrics and further examined through the flux responses to their environmental controls. Our results suggest that CLASSIC tended to overestimate GPP and ER among all stands. Model performance regarding the energy fluxes (i.e., LE and H) varied greatly among the stands and exhibited a moderate correlation with latitude. We identified strong relationships between simulated fluxes and their environmental controls except for H, thus highlighting current strengths and limitations of CLASSIC. 

Background The Drought Code (DC) of the Canadian Fire Weather Index System (CFWIS) has been intuitively regarded by fire managers in Alaska, USA, as poorly representing the moisture content in the forest floor in lowland taiga forests on permafrost soils. Aims The aim of this study was to evaluate the DC using its own framework of water balance as cumulative additions of daily precipitation and substractions of actual evaporation. Methods We used eddy covariance measurements (EC) from three flux towers in Interior Alaska as a benchmark of natural evaporation. Key results The DC water balance model overpredicted drought for all 14 site-years that we analysed. Errors in water balance cumulated to 109 mm by the end of the season, which was 54% of the soil water storage capacity of the DC model. Median daily water balance was 6.3 times lower than that measured by EC. Conclusions About half the error in the model was due to correction of precipitation for canopy effects. The other half was due to dependence of the actual evaporation rate on the proportional ‘fullness’ of soil water storage in the DC model. Implications Fire danger situational awareness is improved by ignoring the DC in the CFWIS for boreal forests occurring on permafrost. 

Abstract The Arctic–Boreal Zone is rapidly warming, impacting its large soil carbon stocks. Here we use a new compilation of terrestrial ecosystem CO₂fluxes, geospatial datasets and random forest models to show that although the Arctic–Boreal Zone was overall an increasing terrestrial CO₂sink from 2001 to 2020 (mean ± standard deviation in net ecosystem exchange, −548 ± 140 Tg C yr⁻¹; trend, −14 Tg C yr⁻¹;P < 0.001), more than 30% of the region was a net CO₂source. Tundra regions may have already started to function on average as CO₂sources, demonstrating a shift in carbon dynamics. When fire emissions are factored in, the increasing Arctic–Boreal Zone sink is no longer statistically significant (budget, −319 ± 140 Tg C yr⁻¹; trend, −9 Tg C yr⁻¹), and the permafrost region becomes CO₂neutral (budget, −24 ± 123 Tg C yr⁻¹; trend, −3 Tg C yr⁻¹), underscoring the importance of fire in this region. 

Abstract We examined the seasonality of photosynthesis in 46 evergreen needleleaf (evergreen needleleaf forests (ENF)) and deciduous broadleaf (deciduous broadleaf forests (DBF)) forests across North America and Eurasia. We quantified the onset and end (Start_GPPand End_GPP) of photosynthesis in spring and autumn based on the response of net ecosystem exchange of CO₂to sunlight. To test the hypothesis that snowmelt is required for photosynthesis to begin, these were compared with end of snowmelt derived from soil temperature. ENF forests achieved 10% of summer photosynthetic capacity ∼3 weeks before end of snowmelt, while DBF forests achieved that capacity ∼4 weeks afterward. DBF forests increased photosynthetic capacity in spring faster (1.95% d⁻¹) than ENF (1.10% d⁻¹), and their active season length (End_GPP–Start_GPP) was ∼50 days shorter. We hypothesized that warming has influenced timing of the photosynthesis season. We found minimal evidence for long‐term change in Start_GPP, End_GPP, or air temperature, but their interannual anomalies were significantly correlated. Warmer weather was associated with earlier Start_GPP(1.3–2.5 days °C⁻¹) or later End_GPP(1.5–1.8 days °C⁻¹, depending on forest type and month). Finally, we tested whether existing phenological models could predict Start_GPPand End_GPP. For ENF forests, air temperature‐ and daylength‐based models provided best predictions for Start_GPP, while a chilling‐degree‐day model was best for End_GPP. The root mean square errors (RMSE) between predicted and observed Start_GPPand End_GPPwere 11.7 and 11.3 days, respectively. For DBF forests, temperature‐ and daylength‐based models yielded the best results (RMSE 6.3 and 10.5 days). 

Environmental observation networks, such as AmeriFlux, are foundational for monitoring ecosystem response to climate change, management practices, and natural disturbances; however, their effectiveness depends on their representativeness for the regions or continents. We proposed an empirical, time series approach to quantify the similarity of ecosystem fluxes across AmeriFlux sites. We extracted the diel and seasonal characteristics (i.e., amplitudes, phases) from carbon dioxide, water vapor, energy, and momentum fluxes, which reflect the effects of climate, plant phenology, and ecophysiology on the observations, and explored the potential aggregations of AmeriFlux sites through hierarchical clustering. While net radiation and temperature showed latitudinal clustering as expected, flux variables revealed a more uneven clustering with many small (number of sites < 5), unique groups and a few large (> 100) to intermediate (15–70) groups, highlighting the significant ecological regulations of ecosystem fluxes. Many identified unique groups were from under-sampled ecoregions and biome types of the International Geosphere-Biosphere Programme (IGBP), with distinct flux dynamics compared to the rest of the network. At the finer spatial scale, local topography, disturbance, management, edaphic, and hydrological regimes further enlarge the difference in flux dynamics within the groups. Nonetheless, our clustering approach is a data-driven method to interpret the AmeriFlux network, informing future cross-site syntheses, upscaling, and model-data benchmarking research. Finally, we highlighted the unique and underrepresented sites in the AmeriFlux network, which were found mainly in Hawaii and Latin America, mountains, and at under- sampled IGBP types (e.g., urban, open water), motivating the incorporation of new/unregistered sites from these groups. 

Abstract Photosynthesis of terrestrial ecosystems in the Arctic-Boreal region is a critical part of the global carbon cycle. Solar-induced chlorophyll Fluorescence (SIF), a promising proxy for photosynthesis with physiological insight, has been used to track gross primary production (GPP) at regional scales. Recent studies have constructed empirical relationships between SIF and eddy covariance-derived GPP as a first step to predicting global GPP. However, high latitudes pose two specific challenges: (a) Unique plant species and land cover types in the Arctic–Boreal region are not included in the generalized SIF-GPP relationship from lower latitudes, and (b) the complex terrain and sub-pixel land cover further complicate the interpretation of the SIF-GPP relationship. In this study, we focused on the Arctic-Boreal vulnerability experiment (ABoVE) domain and evaluated the empirical relationships between SIF for high latitudes from the TROPOspheric Monitoring Instrument (TROPOMI) and a state-of-the-art machine learning GPP product (FluxCom). For the first time, we report the regression slope, linear correlation coefficient, and the goodness of the fit of SIF-GPP relationships for Arctic-Boreal land cover types with extensive spatial coverage. We found several potential issues specific to the Arctic-Boreal region that should be considered: (a) unrealistically high FluxCom GPP due to the presence of snow and water at the subpixel scale; (b) changing biomass distribution and SIF-GPP relationship along elevational gradients, and (c) limited perspective and misrepresentation of heterogeneous land cover across spatial resolutions. Taken together, our results will help improve the estimation of GPP using SIF in terrestrial biosphere models and cope with model-data uncertainties in the Arctic-Boreal region. 

Abstract Tundra and boreal ecosystems encompass the northern circumpolar permafrost region and are experiencing rapid environmental change with important implications for the global carbon (C) budget. We analysed multi-decadal time series containing 302 annual estimates of carbon dioxide (CO₂) flux across 70 permafrost and non-permafrost ecosystems, and 672 estimates of summer CO₂flux across 181 ecosystems. We find an increase in the annual CO₂sink across non-permafrost ecosystems but not permafrost ecosystems, despite similar increases in summer uptake. Thus, recent non-growing-season CO₂losses have substantially impacted the CO₂balance of permafrost ecosystems. Furthermore, analysis of interannual variability reveals warmer summers amplify the C cycle (increase productivity and respiration) at putatively nitrogen-limited sites and at sites less reliant on summer precipitation for water use. Our findings suggest that water and nutrient availability will be important predictors of the C-cycle response of these ecosystems to future warming. 

Abstract. Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr−1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr−1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×103 ± 0.9×103 km3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.